Infrared features of unquenched finite temperature lattice Landau gauge QCD 

Sadataka Furu£] 

School of Science and Engineering, Teikyo University, 320-8551 Japan. 

Hideo NakajimgQ 

Department of Information Science, Utsunomiya University, 321-8585 Japan. 

(Dated: February 1, 2008) 

The color diagonal and color antisymmetric ghost propagators slightly above T c of Nf = 2 MILC 
24 3 x 12 lattices are measured and compared with zero temperature unquenched Nf =2+1 MILC C 
20 3 x 64 and MILC/ 28 3 x 96 lattices and zero temperature quenched 56 4 /3 = 6.4 and 6.45 lattices. 
The expectation value of the color antisymmetric ghost propagator <j> c {q) is zero but its Binder 
cumulant, which is consistent with that of N% — 1 dimensional Gaussian distribution below T c , 
decreases above T c . Although the color diagonal ghost propagator is temperature independent, the 
I 1 norm of the color antisymmetric ghost propagator is temperature dependent. The expectation 
value of the ghost condensate observed at zero temperature unquenched configuration is consistent 
with in T > T c . 

We also measure transverse, magnetic and electric gluon propagator and extract gluon screening 
masses. The running coupling measured from the product of the gluon dressing function and the 
ghost dressing function are almost temperature independent but the effect of A 2 condensate observed 
at zero temperature is consistent with in T > T c . 

The transverse gluon dressing function at low temperature has a peak in the infrared at low 
temperature but it becomes flatter at high temperature. The magnetic gluon propagator at high 
momentum depends on the temperature. These data imply that the magnetic gluon propagator 
and the color antisymmetric ghost propagator are affected by the presence of dynamical quarks and 
there are strong non-perturbative effects through the temperature-dependent color anti-symmetric 
ghost propagator. 

PACS numbers: 12.38.Gc, 12. 38. Aw, ll.10.Gh, 11.15.Ha, 11.15.Tk, 11.30.Rd 



I. INTRODUCTION 

As a condition of the color confinement in the in- 
frared region, Kugo-Ojima criterion[l| and the Gribov- 
Zwanziger scenario|2j, |3| are well-known. In these theo- 
ries, infrared divergence of the ghost propagator is the 
essential ingredient of the emergence of the string-like 
inter-quark potential. In finite temperature SU(2) lattice 
Coulomb gauge simulation, linearly rising color-Coulomb 
potential was observed to remain above the deconfine- 
ment temperature T c [ji|. The color-Coulomb potential 
is defined by the ghost propagator and the color diago- 
nal ghost propagator is found to be essentially temper- 
ature independent. The temperature dependence of the 
color antisymmetric ghost porpagator was not explored. 
Coulomb gauge is a non-covariant gauge and the momen- 
tum in the time direction is not affected by the gauge 
fixing. In finite temperature, the momentum in the time 
direction is interpreted as the Matsubara frequency with 
appropriate boundary conditions. 

The Landau gauge is a covariant gauge, and in the 
analysis of the quenched and unquenched lattice Lan- 
dau gauge simulation [f| we analyzed zero temperature 
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Nf = 2 + 1 unquenched configuration of the SU(3) MILC 
collaboration of 1) 20 3 x 64 (3 imp = 6.76,6.83, which is 
denoted as MILC C , since realtively coarse lattice of spac- 
ing a = 0.12fm is used, and 2) 28 3 x 96 /3 lmp = 7.09, 7.11, 
which is denoted as MILC/, since relatively fine lattice 
of spacing a — 0.09fm is usedQ. 

We observed that the Kugo-Ojima parameters of 
MILCc and MILC/ are consistent with 1 while the 
quenched configuration of 56 4 lattice remained about 
0.8(1, @- The Binder cumulant of the color anti- 
symmetric ghost propagator of the zero temperature 
quenched SU(2) and unquenched Nf =2+1, SU(3) con- 
figurations of the MILC collaboration were consistent 
with those the N% — 1 dimensional Gaussian distribution, 
where N c is the number of colors. The dynamical quark 
has the effect of quenching randomness of the system[l(| ■ 

In this paper, we extend the analysis of the color 
anti-symmetric ghost propagator to the quenched SU(3) 
56 4 lattices of f3 = 6.4@ and 6.45@ and the finite 
temperature Nf = 2 unquenched SU(3) 24 3 x 12 con- 
figurations of j9 = 5.65,5.725,5.85, produced by the 
MILC collaborationfll]. We denote the finite temper- 
ature configurations as MILC/*. The ghost propagator 
contributes in the off-shell quark gluon vertex via Ward- 
Slavnov- Taylor identity [Hj], and so we expect the screen- 
ing masses of the gluons produced by the quark loops and 
ghost loops would be affected by the ghost propagator. 
We study the magnetic and the electric screening mass of 
the gluon of the MILC finite temperature configurations. 
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In the analysis of [ll|, the three configurations corre- 
spond to the temperature T = 143, 172.5 and 185 MeV p , 
respectively where subscript p means that the scale is 
fixed from the mass of the p meson m p = 770MeV. The 
temperature T c at which the chiral susceptibility shows a 
peak, indicating the cross over to the deconfinement was 
assigned to be about 140MeV. Since the standard Wilson 
plaquette action was used in the production of the gauge 
configuration, the flavor symmetry was broken and the 
ratio of the p mass to pion mass is larger than the phys- 
ical value Recent simulation with Asqtad action [14j 
suggests that T c ~ 170MeV, consistent with the value of 
that of Nf = 2 improved Kogut-Susskind (KS) fermion 
T c ~ 173 ± 8MeV [l5|. A systematic comparison of the p 
meson mass for improved staggered actions in quenched 
approximation is given in [16| . Recently, [17| claims that 
the transition temperature depends on which physical 
quantity one measures, and that the T c defined by the 
chiral susce ptib ility is 151(3)MeV and consistent within 
errors with [11| . They showed that the T c defined by the 
strange quark chiral susceptibility and the Polyakov loop 
are about 175MeV and consistent with [l5l |. 

Since the continuum limit of the mass of the vec- 
tor meson would not depend on the temperature near 
T c 18], it would be natural to assign the bare lattice p 
mass by about 20% heavier and shift the temperature 
by the same amount. We leave a more accurate assign- 
ment of the temperature scale of the MILC/t to a fu- 
ture study and since the disconnected chiral susceptibil- 
ity measured by the configurations shows a clear peak 
at around 140MeV p , we assign the (3 = 5.65,5.725 and 
5.85 data, which correspond to T = 143, 172.5 and 185 
MeV p , by T/T c = 1.02,1.23 and 1.32, respectively. We 
compare quenched and unquenched ghost propagator of 
zero temperature and finite temperature and investigate 
the role of quarks on the gluon field and its dependence 
on the temperature. 

The organization of the paper is as follows. In the 
Sect. 2, we show the results of the color diagonal ghost 
propagator of the quenched 56 4 and unquenched finite 
temperature MILC configurations. The corresponding 
color antisymmetric ghost propagators are shown in the 
Sect. 3. The Kugo-Ojima parameters, the transverse, 
magnetic and electric gluon propagator and the QCD 
running coupling of the finite temperature unquenched 
configurations are shown in the Sects. 4, 5 and 6, respec- 
tively. Conclusions and a discussion are given in the 
Sect.7. 



variant derivative, as 

FT[D a G b (x,y)] = FT (tr(A a l {{Mp})- 1 } xy A b )), 

= 5 ab D G {q 2 ). (1) 

where U is the link variable obtained by the Landau 
gauge fixing. In all the simulation in this work we adopt 
the \ogU definition of the gauge field, and the SU(3) 
color matrix A is normalized as trA a A b = 5 ab . Number 
of samples is about 10 each in the 56 4 lattices and about 
100 each in the 24 3 x 12 lattices. 

In the usual Monte Carlo simulation, it is necessary to 
make an average over many samples, but the color diago- 
nal ghost propagator and the color anti-symmetric ghost 
propagator of unquenched configurations are almost in- 
dependent of samples. Since we adopt the cylinder cut, 
i.e. select the momentum along the diagonal direction in 
momentum space, and take into account translation in- 
variance and rotational symmetry of the lattice, the error 
bar of the ghost propagator of one sample is already not 
so large, and an order of 10 samples is enough to get the 
expectation values, although one should be cautious to 
the appearance of exceptional samples Q. 

The color diagonal ghost propagator is defined as 

° G{q) = Nf—lV 

xtr(S ab ((A a cosqx\f b (x)) + (A a smqx\f b (x)))) 

= G(g 2 )/g 2 



p b (x) with p b (x) = — — A b cosga; and — —A b sinqx, re- 

Vv Vv 



Here fc(x) and f s h {x) are the solution of Mf b (x) — 
p b (x) with 
spcctively. 

A. Quenched SU(3) 56 4 lattice 

The color diagonal ghost propagators of 56 4 quenched 
SU(3) are shown in [6j, but we show the data for a com- 
parison with the color anti-symmetric ghost propagator 
in the next subsection. The FIG. |T] is the ghost dressing 
function, of the 56 4 j3 = 6.45 configurations produced 
by the Monte Carlo simulation and subsequently Landau 
gauge fixed. 

The scale of the j3 — 6.4 and 6.45 configurations are 
fixed as in TABLE [fl using the formula in [19]. The 
ghost propagator is infrared divergent, but its singularity 
is weaker than q~ A . 



II. THE COLOR DIAGONAL GHOST 
PROPAGATOR 

The ghost propagator is defined by the Fourier trans- 
form(FT) of the expectation value of the inverse Faddeev- 
Popov(FP) operator M. = —d^D^, where D p is the co- 



B. MILC finite temperature 24 3 x 12 lattice 

The color diagonal ghost dressing function of MILC 
finite temperature configurations are shown in FIG. [2] 
We observe that the three data of different temperatures 
scale as shown in the figure when each lattice spacing a 
is chosen as in the ref. [ill . 
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FIG. 1: The ghost dressing function of quenched 56 4 config- 
urations of P = 6.45 as the function of log 10 [q(GeV)]. 



FIG. 2: The ghost dressing function of MILC/t configurations 
of T/Tc = 1.02(blue diamonds), T/T c = 1.23(red stars) and 
T/T c — 1.32(green triangles) as the function of log 10 [q(GeV)]. 



TABLE I: The (3 = 2N c /g 2 {(3 irnp = 5/3 x (3 for MILC/ 
and MILC C ), the inverse lattice spacing l/o, lattice size 
and lattice length(fm) of configurations investigated in this 
paper. Subscripts c and / of MILC correspond to coarse 
lattice(a=0.12fm) and fine lattice(a=0.09fm). The MILC/ t 
corresponds to the Nf = 2 finite temperature configurations 
and we assign /3 = 5.65, 5.725 and 5.85 data to T/T c — 
1.02, 1.23 and 1.32, respectively. 





f3imp / f3 


am u d/ arris Nf 


l/a(GeV) L B L t 


oLs (fm) 


quench 


6.4 





3.66 


56 56 


2.96 




6.45 





3.87 


56 56 


2.94 


MILC C 


6.83 


0.040/0.050 2+1 


1.64 


20 64 


2.41 




6.76 


0.007/0.050 2+1 


1.64 


20 64 


2.41 


MILC/ 


7.11 


0.0124/0.031 2+1 


2.19 


28 96 


2.52 




7.09 


0.0062/0.031 2+1 


2.19 


28 96 


2.52 


MILC /t 


5.65 


0.008 2 


1.716 


24 12 


2.76 




5.725 


0.008 2 


1.914 


24 12 


2.47 




5.85 


0.008 2 


2.244 


24 12 


2.11 



III. THE COLOR ANTISYMMETRIC GHOST 
PROPAGATOR 

The color anti-symmetric ghost propagator is defined 



: (<z) 



i i 

Jfv 

cabc 1 1 a a 



xtr(r bc ((A a co Sq x\f b s (x)) - (A a smqx\f b c (x)))) 

where the outer-most bracket means the ensemble aver- 
age. 

In a theory based on the Curci-Ferrari gauge and in 
its extension to the Landau gauge, a parameterization 
of the color anti-symmetric ghost propagator with use of 
the ghost condensate parameter v was proposed [20l|. 

A simulation in the SU(2) lattice Landau gauge using 
a parameterization 
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N? - 1 ^ 



^\r( q )\- r/L 



9 4 + v l 



(2) 



is performed in [21(. In [2(| the ghost propagator in the 
condensed vacuum in SU(2) is expressed as 



a b\ 



2 sab 



(C C 



(3) 



The term r/L 2 in eq.([2]) is a term introduced to incorpo- 
rate the finite lattice size effect. We applied the method 
to zero temperature unquenched SU(3) configurations of 

MILCc and MILC/0,0]- In SI, nnite size effects ex- 
pressed by r/L 2 is disentangled from the fitting of the lat- 
tice data of color antisymmetric ghost propagator in high 
momentum region, where the perturbative QCD (pQCD) 
is a good approximation, as 



L 2 



where the denominator in 



cos(7rg/L) 

I0 a (9)l 



\r(q)\ 



(4) 



is the factor that 



cos(7rq/L) 

appears in the vertex function of the lattice perturba- 
tion theory and q takes integer values 0,1, ■■■L/2. In 
the present work, since the high momentum data points 
where the enhancement due to 1/ cos factor is significant 
are few, i.e. the enhancement of the highest momentum 
point of FIGH] is 0.64, that of the next to highest mo- 
mentum point is 0.35 and the rest are less than 0.2 in 
the log 10 scale, we fix the parameter r and v by the fit 
of |0 a (g)| using eq.© including the whole momentum 
region except the highest momentum point. 

This treatment of the finite size effect on the lattice 
data is not without ambiguity, but our aim is to search 
qualitative differences as the temperature or the masses 
of quarks are changed. 

Another quantity that characterizes the system is the 
Binder cumulant of the color antisymmetric ghost prop- 
agator defined as 



U(q) = 1 



Wg) 4 ) 

mi) 2 ) 2 ' 



A simulation of SU(2) lattice Landau gauge obtained 
U ~ 0.44, almost independent of the momentum. This 
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FIG. 3: log 10 \<j>(q)\ as the function of log 10 g(GeV) of FIG. 4: log 10 \<j>(q)\ as the function of log 10 g(GeV) of MILC/ 



quenched (3 = 6.45, 56 lattice. 



and the fit using r = 134 and v = 0.026GeW 



value is compatible with that of the three dimensional 
Gaussian distribution (Toj and the analysis of SU(3) 
MILC C and MILC/ showed that U is compatible with 
that of eight dimensional Gaussian distribution. We ex- 
tend these analyses to large quenched lattices and the 
finite temperature configurations. 



A. Quenched SU(3) 56 4 lattice 

The absolute value of the color antisymmetric ghost 
propagator of quenched 56 4 lattice in the infrared region 
is about 3 orders of magnitude smaller than that of the 
color diagonal ghost propagator and the values are sam- 
ple dependent. Results of (3 = 6.45, 56 4 10 samples are 
shown in FIG [3] 

Due to this sample dependence, quite different from 
the case of unquenched 20 3 x 64 lattices [10| and 28 3 x 
96@, the Binder cumulant of the color antisymmetric 
ghost propagator of quenched configurations is noisy due 
to large (4>(q) 4 ) as compared to ((j>(q) 2 ) 2 . The random- 
ness of the color anti-symmetric ghost propagator of 
quenched configurations is large and the Binder cumu- 
lant becomes unstable. 



B. MILC zero temperature 28 x 96 lattice 

The color antisymmetric ghost propagator of MILC/ 
is shown in FIG. |U The fitting parameters are given in 
1- 

The expectation value of v is small but the condition 
that the fitted curve passes the lowest momentum point 
of FIGHwithin the error bar gives v = 0.026(6)GeV 2 . It 
suggests that the presence of BRST partner of A 2 con- 
densate at zero temperature. 

The Binder cumulant of the zero temperature Nf = 
2 + 1 MILC C lattice is reported in [loj] and that of MILC / 
is reported in Q. We observed that the mass function 
of m = 27.2MeV quark propagator of [3i mp — 7.09 
with bare mass combination mp = 27.2Me V/68MeV 



shows an anomalous behavior in q < lGeV region, and 
that in the same region Binder cumulant U(q) shows 
an anomalous behavior, although the mass function of 
rno = 68MeV does not show the anomaly. The non-QCD 
like behavior of staggered quarks calculated with large 
lattice spacing a and small bare mass too is reported 
in [23]. Since no anomaly was observed in = 7.09 
too = 13.6MeV/68MeV @,E3|, effects of the relative size 
of the s— quark mass v. ud— quark mass and the number 
of Nf are suggested. 

Thus we extend the analysis to MILC configuration 
7.18, 28 3 x 96 lattice with the bare 



of Nf = 3, I3 imp 
quark mass too = 

20 3 x 64 lattice with the bare quark mass toq = 0.02 [31 1 



0.031 and that of Nf = 2, # 



7.20, 



A result of the Binder cumulants of the Nf — 3 and 
the Nf = 2 (50 samples) are 0.66(1), i.e. {${q) A ) is 
close to ((f>(q) 2 ) 2 . When the bare masses of the quarks 
are the same, the system possesses the self averaging 
property [32|], and when they are different as in Nf = 2+1, 
the system lacks this property. The parameter v fitted 
from Nf = 2 (50 samples) is found to be consistent with 
0. 



C. MILC finite temperature 24 3 x 12 lattice 

The fluctuation of the color antisymmetric ghost 
propagator around the expectation value was almost 
Gaussian in the case of zero temperature unquenched 
configurations [9I fioj] . The logarithm of the absolute value 
of the color antisymmetric ghost propagator of finite tem- 
perature unquenched configurations as a function of the 
logarithm of the momentum is shown in the FIG. [5] In 
contrast to the color diagonal ghost propagator, the ab- 
solute value of the color antisymmetric ghost propagator 
depends on the temperature. The temperature depen- 
dence of the scale of cj) defined by r of eq. ((4]) can be ex- 
pressed as roughly r ~ 5.49a~ 4 23 , where a is the lattice 
spacing at each temperature in the unit of GeV. 

The fitting parameters of the color anti-symmetric 
ghost propagator are given in TABLEQl] The condensate 
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FIG. 5: log 10 \</>(q)\ as the function of log 10 q(GeV) of MILC /t 
of T/Tc = 1.02(blue diamonds), T/T c = 1.23(red stars) and 
T/T c = 1.32 (green triangles). 



FIG. 6: The Binder cumulant of the color antisymmetric 
ghost propagator of MILC Nf = 2 configurations of T/T c = 
1.23(red stars). 



parameter v of finite temperature is difficult to assign 
since the lattice spacing is relatively large and difficult 
to detect the infrared bending behavior in the color an- 
tisymmetric ghost propagator, but they are consistent 
with 0. 

TABLE II: The fitted parameters r, z and v of the color an- 
tisymmetric ghost propagator \(/>(q)\ of MILC C , MILC/ and 
MILC finite temperature. Two values of U of MILC/ corre- 
spond to the average below q — lGeV and the average above 
lGeV, respectively. 



Pimp/ 'P mo(MeV) 


r 


z v U 


6.76 


11.5/82.2 


37.5 


3.90 0.02(1) 0.53(5) 


6.83 


65.7/82.2 


38.7 


3.85 0.01(1) 0.57(4) 


7.09 


13.6/68.0 


134 


3.83 0.026(6) 0.57(4)/0.56(l) 


7.11 


27.2/68.0 


112 


3.81 0.028(8) 0.58(2)/0.52(l) 


5.65 


12.3 


54.4 


4.01 0.0 0.580(13) 


5.725 


12.8 


88.3 


3.95 0.0 0.571(4) 


5.85 


15.0 


165.9 


3.93 0.0 0.558(2) 



The Binder cumulant of (f> is almost independent of 
the momentum except the lowest momentum point. We 
show a typical example of T/T c = 1.23 in FIG. M 

The temperature dependence of the average of U(q) 
excluding the lowest momentum point is shown in FIG. [7] 
The average U (T) decreases monotonically as a function 
of T from that of the Gaussian distribution (dashed line) 
at high temperature. 



IV. THE KUGO-OJIMA COLOR 
CONFINEMENT PARAMETER 

The Kugo-Ojima parameter is defined by the two point 
function of the covariant derivative of the ghost and the 
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FIG. 7: Averages over momenta excluding the lowest momen- 
tum point of the Binder cumulants of MILC finite tempera- 
ture configurations. T/T c — 1.02(blue diamonds), T/T c = 
1.23(red stars) and T/T c = 1.32(green triangles). 



commutator of the antighost and gauge field 




(5) 

Kugo and Ojima[l| showed that u(0) = — 1 is a condi- 
tion of the color confinement. ZwanzigerQ defined the 
horizon function h that is related to the Kugo-Ojima pa- 
rameter c = —u(Q) as follows. 

e~ tq{x ~ y) ( tr (V^-^- (-£„) A 6 ) \ 

x,y \ > Sxy/ 

= G, u (q)5 ab = (I) ^s ab - (V - u ab , 

where, with use of the covariant derivative D^fJJ) 
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4>{x + n) - <j>{x), 
_ adjA x ^/2 

t&nh(adjA X41 /2) 
Using the definition 



(x + ft) + <j){x) 



and 



S{U XI1 ) 



= (^tr(A°t5(C7^)A Q ))/{(iV c 2 -l)y}, 



the horizon condition reads lim G MM (g) — e = 0, and the 



left hand side of the condition is 



(|) + («-l)< 



[d — l)h where h = c — — and dimension d = 4, and it 

follows that h = — > horizon condition, and thus the 
horizon condition coincides with Kugo-Ojima criterion 
provided the covariant derivative approaches the naive 
continuum limit, i.e., e/d = 1. 

The Kugo-Ojima parameter is defined by a scalar func- 
tion at vanishing momentum in the continuum theory, 
but in the lattice simulation, we measured the magni- 
tude of the right hand side of the eq.© with \i = v 
polarizations of A v and D^, and then observed in the 
case of asymmetric lattices there is a strong positive cor- 
relation between the magnitude and the lattice size of the 
axis whose directions are perpendicular to the polariza- 
tion. This general features are manifested in the data of 
MILC/ and MILC C , whose time (4th) axis is longer than 
the spatial axes and in the data of MILC/* whose time 
axis is shorter than the spatial axes. 



A. Quenched SU(3) 56 4 lattice 
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FIG. 8: Kugo-Ojima parameter u(0) of quenched 56 4 con- 
figurations of j3 = 6.45(red stars). The data of an excep- 
tional sample is indicated by the green diamond. Polariza- 
tions 1,2,3,4 correspond to x,y,z and t. 
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FIG. 9: Kugo-Ojima parameter w(0) of MILC/ Nf =2 + 1 
KS fermion unquenched configurations of /3i mp = 7. II (green 
diamonds), f3im P = 7.09(red stars) . 



The Kugo-Ojima confinement parameter of quenched 
configuration saturates at about 80% of the expected 
value c = 1. There appear exceptional samples with an 
average consistent with c = 1 within errors. The po- 
larization dependence of the Kugo-Ojima parameter of 
/3 = 6.45 samples shown in FIG. [5] is due to the lack of 
rotational invariance of our system. 

B. MILC finite temperature Nf = 2, 24 3 x 12 lattice 



large when there is a long axis perpendicular to the po- 
larization. The FIG. [TOl shows the temperature and the 
polarization dependence of the Kugo-Ojima parameter. 

The quenched configurations and high temperature 
configurations show larger deviation from u(0) = — 1, 
which would be due to higher randomness of these sys- 
tems. The origin of the randomness in the latter would be 
the thermal fluctuation and that of the former would be 
lack of fermions that quench randomness of the system. 



Shown in TABLE IIIII are the Kugo-Ojima parame- 
ters and the horizon function deviation of quenched 56 4 
configurations, zero-temperature unquenched MILC C and 
MILC/ configurations and unquenched finite tempera- 
ture MILC ft configurations. 

The Kugo-Ojima parameters of unquenched zero- 
temperature configurations are consistent with the 
theory [8J, |{J, [ltj, and those of unquenched finite temper- 
ature configurations show temperature dependence. 

The FIG. [9] show the dependence on the polarization 
of the Kugo-Ojima parameter. The polarization depen- 
dence of the parameter c is a result of an integration over 
the axes perpendicular to the polarization. It becomes 



V. THE FINITE TEMPERATURE GLUON 
PROPAGATOR 

The fluctuation of the gluon propagators of lattice Lan- 
dau gauge is larger than that of the ghost propagator. 
We measure the gluon propagators of the three temper- 
atures of MILC ft by using about 100 samples each and 
choosing the momenta along the three dimensional (3-d) 
coordinate axes and along the diagonal in the 3-d space. 
The dependence on the choice of the direction is not sig- 
nificant except a slight suppression of the q = (1, 1, 1, 0), 
as compared to that with the same magnitude but along 
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quench 


6.4 




0.827(27) 


0.954(1) 


-0.12 




6.45 




0.814(89) 


0.954(1) 


-0.14 


MILC C 


6.76 


1.04(11) 0.74(3) 


0.97(16) 


0.9325(1) 


0.03(16) 




6.83 


0.99(14) 0.75(3) 


0.93(16) 


0.9339(1) 


-0.00(16) 


MILC/ 


7.09 


1.06(13) 0.76(3) 


0.99(17) 


0.9409(1) 


0.04(17) 




7.11 


1.05(13) 0.76(3) 


0.98(17) 


0.9412(1) 


0.04(17) 


MILC/t 


5.65 


0.72(13) 1.04(23) 


0.80(21) 


0.9400(7) 


-0.14(21) 




5.725 


0.68(15) 0.77(16) 


0.70(15) 


0.9430(2) 


-0.24(15) 




5.85 


0.63(19) 0.60(12) 


0.62(17) 


0.9465(2) 


-0.33(17) 



TABLE III: The Kugo-Ojima parameter of the quenched 56 4 lattice and that of the unquenched KS fermion 
(MILCc, MILC/, MILC/ t ). c x is the polarization along the spatial directions, c t is that along the time direction, c is the 
weighted average of c x and Ct i.e. (3c x + ct)/4, e/d is the trace divided by the dimension and h is the horizon function 
deviation. 





-0.2 
-0 . 4 
-0 . 6 
-0 . 8 
-1 
-1.2 
-1 . 4 



\ I \\ 



12 3 
polarization 



FIG. 10: Kugo-Ojima parameter u(0) of MILC/t configura- 
tions of T/T c = 1.02(blue diamonds), T/T c = 1.23(red stars) 
and T/T c = 1.32(green triangles). 



the coordinate axes. We fix the normalization of the 
gluon propagator such that the running coupling denned 
by the product of the gluon dressing function and the 
ghost dressing function squared as defined in Sect I VII in 
the high momentum region agrees with the pQCD result 
and that the rescaling factor for the ghost dressing func- 
tion and the gluon dressing function are the same. The 
renormalization factors defined by the highest momen- 
tum point along the coordinate axes of the three temper- 
atures are 0.48(1). 

In Fig[Tl]the transverse gluon dressing functions of the 
three temperatures taken at momenta along the diagonal 
direction in the 3-d momentum space after the rescaling 
are shown. The asymptotic value is about 2, since we 
normalized TrA a A b — 5 ab . Larger infrared enhancement 
of the transverse gluon dressing function for low temper- 
ature agrees with a recent result of quenched SU(2) finite 
temperature simulation [23j. 

The temperature dependence of unquenched trans- 
verse gluon propagator and that of the magnetic gluon 
propagator G m (z), and the temperature dependence of 
the quenched SU(2) transverse gluon propagator [23j is 
similar to that of unquenched SU(3) electric gluon prop- 



agator G e (z), which is shown in Fig[T4l and [T3l respec- 
tively. 

The transverse gluon propagator at the zero momen- 
tum is finite as shown in Fig|T2j The value of Da(0) 
of MILC/, (3 = 7.09, 28 3 x 96 is about 15 GeV- 2 (|. 
The -0^(0) decreases monotonically as T decreases to 
0. Whether it vanishes or remains finite is an important 
problem for fixing the nature of the infrared fixed point 
of the running coupling. There is an argument that the 
infrared non-vanishing of the Landau gauge gluon prop- 
agator is a finite size effect. Our lattice volume of low 
temperature data is larger than the high temperature 
data, but the running coupling which is calculated by 
the gluon propagator and shown in Sect. VI excludes 
large temperature dependence in the finite size effects. 
A correction of the global scaling factor is, however, not 
excluded. 




2 3 
q(GeV) 



FIG. 11: The transverse gluon dressing function Z(q) of 
MILC/t Nf = 2 KS fermion unquenched configuration of 
T/T c = 1.02(blue diamonds), T/T c = 1.23(red stars) and 
T/T c = 1.32(green triangles). 



We observe that the fluctuation of the gluon propaga- 
tor at the zero momentum increases as the temperature 
increases.® 

The one-loop off-shell contribution to the quark-gluon 
vertex, is related to the quark-quark-ghost-ghost ampli- 
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1 1.1 1.2 1.3 2 4 6 8 10 12 

T/T c z 



FIG. 12: The transverse gluon propagator of MILC/t Nf = 2 
KS fermion unquenched configuration of T/T c = 1.02(blue 
diamonds), T/T c = 1.23(red stars) and T/T c = 1.32 (green 
triangles) at q = in the unit of GeV -2 . 



FIG. 13: Electric gluon propagator of MILC /t TV/ = 2 KS 
fermion unquenched configuration of T/T c = 1.02(blue dia- 
monds), T/T c = 1.23(red stars) and T/T c = 1.32(green trian- 
gles). 



tude via Ward-Slavnov- Taylor identity [T^|. The screen- 
ing mass of the gluon would be affected by the ghost 
propagator. The correlation function of the gauge fields 
in the position space is defined as [24[ 

D Afi (z) = (A,(z)A^0)) } 

where A^z) = T, X0 , X1<X2 A„(x , £). 

The screening mass of the gluon is produced by the 
quark-, gluon- and ghost-loops. The electric screening 
mass m e is defined by 

G e (z) = Dao(z) ~ exp[-m e z / N t T] 

The magnetic screening mass m m is defined by 

G m (z) = (D A1 (z) + D A2 (z))/2 ~ exp[-m m z/N t T] 

We fit the propagator G e (z) and G m (z) by 

D A (z) = Acosh[(m/N t T){z - N x /2)}. 

The mass depends on the region of z used to fit the data 
and Karsch et al. proposed to fit data of z near N x /A. [25| 
We fit the data from z — 5 to 9. When we fit the data 
from z — to z = 6, the masses become smaller by about 
30-40%. Since deviations from the cosh functional form 
become larger in this region, we adopt the fit from z = 5 
to 9. 

TABLE IV: The electric and the magnetic screening mass of 
MILC/ t . 

13 T/T c m e /T m m /T 
5.65 1.02 3.42(27) 3.48(48) 
5.725 1.23 3.22(40) 2.90(20) 
5.85 1.32 3.14(33) 2.31(22) 

In high temperature, pQCD suggests that 
^(x.g(T),^W(T) 








2 4 6 8 10 12 



z 



FIG. 14: Magnetic gluon propagator of MILC/t N f = 2 KS 
fermion unquenched configurations of T/T c = 1.02(blue dia- 
monds), T/T c = 1.23(red stars) and T/T c = 1.32(green trian- 
gles). 



where g 2 (T)/47r is the running coupling at temperature 
T and momentum. Our electric screening mass and the 
magnetic screening mass are close to those of quenched 
SU(2) TO e /T = 2.484(52) measured by 0. The mag- 
netic screening mass in the case of SU(2) using the data 
in the region of T > 2T C was m m /T = 0.466(25)g 2 (T) in 
the two-loop pQCD calculation with A m = 0.262(18)T C . 
In the case of quenched SU(3), obtained using the 
data in the region of T > 1.5T C , m e /T = 1.69(4) ff (T) 
and m m /T = 0.549(16)g 2 (T). Since the electric and the 
magnetic gluon propagator of [26| are normalized to 1 
at z — and the critical temperature of the quenched 
configuration T c ~ 269 ± lMeV[l]| is much higher than 
that of the unquenched configuration, we cannot com- 
pare quantitatively their data with ours, but their data 
oiT/T c — 1.32 are consistent with ours within errors. We 
do not observe suppression of m e /T near T c . Whether 
the discrepancy is due to the presence of dynamical 
quarks is left for the future study. Discrepancy of about 
factor 6 in the m m /T of SU(3) near T c from the extrap- 
olation of the pQCD results in T > 1.5T C region implies 
breakdown of the perturbation series near T c . 
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FIG. 15: The magnetic screening mass m m /T (stars) and the 
electric screening mass m e /T of MILC/t Nf = 2 KS fermion 
unquenched configurations (triangles) . 



VI. THE FINITE TEMPERATURE QCD 
RUNNING COUPLING 

In [3(| the calculation of the running coupling in the 
Dyson-Schwinger equation(DSE) at zero temperature [33| 
was extended to below T c and the running coupling for- 
go = 2nnT is written as 

a(n,p 2 ,T) = a( f i 2 )G 2 (n,p 2 ^ 2 ,T)Z(n,p 2 ^ 2 1 T) (6) 

where G(n,p 2 , /j? ,T) is the ghost dressing function and 
Z(n,p 2 , /j 2 , T) is the gluon dressing function. Apply- 
ing this method to T ~ T c , we measure the QCD run- 
ning coupling in the MOM scheme as the product of 
the gluon dressing function and the ghost dressing func- 
tion squared, as a function of the momentum q along 
the spatial axes. We set n equals in the analysis. 
We normalize the magnitude to the result of four-loop 
pQCD[l, H3, [H, [H[ at the highest momentum point of 
/3 = 5.725. A DSE result of zero-temperature quenched 
configuration with the infrared exponent of the ghost 
propagator n = 0.5 (30| . In [fjj we compared the run- 
ning coupling in the quenched lattice simulation with the 
parametrization suggested by the DSE analysis as 



a s {q 2 



1 



2/A2 m 0000 • # 



c + (q 2 /A 2 QCD ) 
1 



4tt 



I 2 tfl 



log(g 2 /A 2 
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) (q 2 /A 2 



QCD 
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'A 2 

^QCD 
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o _ 4tt 8?rft log(i) 
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(8) 



where ft = (102 - 
- ,2857 5033 A7 

P2 = (— — ^/ 



fJV»/2, 



and 



ft = ( 



2 

149753 
6 



18 



325 

1a 



/V 2 )/4, 



+ 3564C(3) + (- 



1078361 6508 



,50065 6472 



C(3))/V 2 



1093 



162 
Nf)/8. 



27 



C(3))JV> 



162 81 

Here the scale in MOM scheme is defined as 13 



A = A 



= (70/3-22JV / /9)/22 



and A^ is 0.259(22) GeV for N f = 2 and 0.252(10) 
GeV for TV/ = were obtained by fitting the data in 
the continuum window 1.8GeV < q <2.3GeV (0.26 < 
log 10 [g(GeV)]<0.37)|El2l- _ 

The pQCD running coupling in the MOM scheme us- 
ing Ajjg = 0.259 GeV is shown by the dotted line. 

Comparing with the result of Nf = 2 + 1 zero- 
temperature [8j, the deviation from the pQCD is smaller. 
Although the presence of ghost condensate at zero tem- 
perature is not excluded, there is no sign of finite v and 
yl 2 at finite temperature. 



4 

3 . 5 

3 

2 . 5 

2 



1.5 

1 

. 5 



4j 



-0.4-0.2 0.20.40 

Log 10 [q(GeV) ] 



6 0.8 



where ft = 11 - fNf. 

The parameters adopted in the fit was Aqcd = 
330MeV, co = 30. The Nf dependence of the running 
coupling is not so large and to make comparison with the 
quenched data, the dashed line used to fit the quenched 
data of Nf = @ is shown in Fig [TBI 

In the pQCD, an approximate inversion of the four- 
loop formula yields the running coupling as a function of 
t = log(g 2 /A 2 ) as follows [13, 0|. 



FIG. 16: The running coupling of the MILC/ t of /3 = 
5.65(blue diamonds), 5.725(red stars) and 5.85(green trian- 
gles) as the function of logq(GeV). Dashed line is the DSE 
result used to fit the quenched zero temperature lattice data 
[(| with k — 0.5 and the dotted line is the pQCD result with 
Nf = 2 of zero temperature. 

The running couplings of three temperatures near T c 
as a function of the momentum q lies roughly on a single 
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curve and there is not strong temperature dependence. 
The relatively large T dependence of the m m /T suggests 
that the non-perturbative effects on the magnetic screen- 
ing mass is important. 

The running coupling of zero temperature Nf = 2 
MILC 20 3 x 64 configuration is roughly the same as that 
of N f = 2 + 1 MILC C presented in Fig. 15 of 8]. In the 
—0.9 < log 10 [q(GeV)] < —0.3 region it decreases mono- 
tonically. We think that the suppression of the running 
coupling at the infrared is due to the singularity of the 
color antisymmetric ghost dressing fuction there, which 
suppresses the color diagonal ghost dressing function, and 
that the suppression is a kind of an artefact. It may im- 
ply a question on the structure of the infrared fixed point 
and on the extensibility from the perturbative region to 
the non-perturbative region of the method of defining the 
running coupling from the product of the color diagonal 
ghost dressing function squared and the gluon dressing 
function. 



VII. DISCUSSION AND CONCLUSION 

We measured the color diagonal and the color antisym- 
metric ghost propagator of quenched, unquenched zero 
temperature and unquenched finite temperature config- 
urations. The quark has the effect of quenching random- 
ness of the ghost propagator and enhancing the trans- 
verse and magnetic gluon propagator above T c . The 
Binder cumulant of the quenched SU(3) color antisym- 
metric ghost propagator is noisy, but that of the un- 
quenched SU(3) color antisymmetric ghost propagator 
is almost constant and independent of the momentum 
except the lowest momentum point. The color diagonal 
ghost propagator of finite temperatures above T c are al- 
most temperature independent, but the scale of the color 
antisymmetric ghost propagator becomes larger as the 
temperature becomes higher. 

We observed that the Binder cumulants of zero tem- 
perature unquenched Nf = 2 + 1 SU(iV c ) configurations 
are consistent with those of N% — 1 dimensional Gaussian 
distributions. The Binder cumulants of finite tempera- 



ture unquenched Nf = 2 + 1 SU(3) configurations devi- 
ate from Gaussian distribution as the temperature rises, 
which may be interpreted as the quark effect of quench- 
ing randomness is reduced in high temperature. The 
anomaly of the KS fermion propagator and the anomaly 
of the momentum dependence of the Binder cumulant 
seems to be correlated. It may indicate that the QCD- 
like region of Nf = 2 + 1 KS fermion system is more 
complicated than that of the Nf = 2 case given in [12] ■ 
The Binder cumulant of unquenched Nf = 2 and 3 which 
are consistent with 0.66 indicates that the ratio of the 
ud— quark mass and the s— quark mass is an important 
ingredient. 

The parameter v introduced to investigate the ghost 
condensate is found to be small and its value depends on 
the temperature and the bare mass of the KS fermion. 

We observed strong non-perturbative effects in the 
magnetic screening mass of the gluon near T c . Since near 
T c , systematic perturbative QCD calculation is impossi- 
ble, it is important to formulate the lattice theory. In 
the quenched finite temperature Landau gauge DSE ap- 
proach [36l . |37| the infrared exponent k and the running 
coupling in the MOM scheme were found to be essentially 
temperature independent below T c . 
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